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SYSTEM AND METHOD OF MEASURING AND MAPPING 
THREE DIMENSIONAL STRUCTURES 



CROSS-REFERENCES TO RELATED APPLICATIONS 

[0001] This application is a continuation-in-part of co-pending U.S. patent application 
10/419,072, filed on 21 April 2003 which is in turn a continuation of U.S. patent application 
09/692,483, filed on 20 October 2000 and which issued on 22 April 2003 as U.S. Patent 
6,550,917, which application in turn claims priority under 35 U.S.C. § 1 19 from U.S. 
Provisional Patent Application No. 60/182,088 filed on 11 February 2000, the entire contents 
of each of which applications are hereby incorporated by reference in their entirety for all 
purposes as if fully set forth herein. This application is also a continuation-in-part of co- 
pending U.S. patent application 10/369,513, filed on 21 February 2003, the entire contents of 
which are hereby incorporated by reference in their entirety for all purposes as if fully set 
forth herein. 



BACKGROUND AND SUMMARY 



[0002] Field. 

[0003] This invention pertains to the field of measurement of surfaces of three-dimensional 
objects, and more particularly, to a system and method of measuring and mapping three- 
dimensional surfaces of spherical and aspherical objects, and of characterizing optically 
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transmissive devices. 
[0004] Description. 

[0005] There are a number of applications where the ability to provide a true measurement 
and map of a three-dimensional surface would be beneficial. There are also a number of 
applications where the ability to precisely and accurately determine characteristics 
(particularly one or more optical properties) of an optically transmissive device would be 
beneficial. 

[0006] For example, contact lenses for correcting vision are produced using contact lens 
molds. Such contact lens molds may suffer not only from surface irregularities (small-scale 
features), but also may include large-area defects. Therefore, it would be desirable to provide 
a system and method of measuring and mapping a three-dimensional surface of a contact lens 
mold. 

[0007] U.S. Patent 6,550,917 ("the '917 patent"), from which this application claims 
priority, describes an improved system and method of compiling a topographic mapping of 
refractive errors in the eye. A disclosed system employs a wavefront sensor, such as a Shack- 
Hartmann wavefront sensor, and includes a combination of an adjustable telescope and a 
range-limiting aperture (RLA) to provide a high dynamic range measurement. 
[0008] However, such a system cannot be directly applied to aspherical measurements of a 
contact lens mold, for example, because in most cases these elements are highly curved 
surfaces that cannot be probed with an injected beam system. Also, for measurements of the 
eye, the aberrations are known to vary according to specific classifications, such as defocus, 
astigmatism, coma or other higher order effects, and thus the ocular instrument could be 
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designed a priori to deal individually with these different, known, types of aberrations. 
[0009] The use of a wavefront sensor for optical metrology of optics, surfaces, optical 
systems, and the characterization of light beams is a well established art. Methods have been 
developed for: characterization of laser beams (e.g., U.S. Patent Number 5,936,720; D. R. 
Neal et al., "Amplitude and phase beam characterization using a two-dimensional wavefront 
sensor," S.P.I.E. 2870 (July 1996) pp.72-82); measurement of optics (e.g., D. R. Neal et al., 
"Wavefront sensors for control and process monitoring in optics manufacture," S.P.I.E. 
2993 , pp. 21 1-220 (1997); J. Pfund et al., "Nonnull testing of rotationally symmetric 
aspheres: a systematic error assessment," Appl. Optics 40(4) pp. 439^446 (1 Feb. 2001); and 
M. Abitol, "Apparatus for mapping optical elements," U.S. Patent 5,825,476); measurement 
of the eye and other ophthalmic systems (U.S. Patent 6,550,917 (the c 917 patent")); and for 
many other applications. If the light incident on a wavefront sensor is the result of reflection 
or scattering light from a surface, then the metrology or other characteristics of the surface 
can be determined. 

[00010] In this art several different technologies have been developed for sensing the 
arriving wavefront of the incident light. Among these technologies are interferometry, Shack- 
Hartmann wavefront sensing, Moire' deflectometry, Shearing interferometry, phase-diversity, 
curvature and other sensors. Each of these different types of sensors has a specific dynamic 
range and accuracy, along with other requirements for proper functioning of the appropriate 
systems. The design and selection of the appropriate phase sensitive instrument depends 
upon the desired characteristics of the measurement system including the desired dynamic 
range and accuracy. 
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[00011] However, in many cases the requirements for dynamic range and/or accuracy exceed 
the ability of a particular measurement technology. For a case where a large dynamic range is 
needed, the accuracy of the measurement may be reduced. If the instrument is designed to 
meet a specific accuracy, then it will often have reduced dynamic range. Thus one or both 
elements of the metrology system (dynamic range or accuracy) must be compromised in order 
to achieve the desired result. For example, Pfund describes a non-null test of an asphere 
using a Shack-Hartmann wavefront sensor approach. However, he carefully limits the 
application to rotationally symmetric systems so that the instrument can be operated within 
the dynamic range of the system. Thus his device is of limited applicability to general 
aspherical shapes. Abitol attempts to overcome this difficulty by marking one lenslet 
differently from the others. This is done by eliminating the central lenslet from the pattern of 
lenslets. Greater dynamic range can be achieved by extrapolating the position of this spot 
outwards from this marked location. However, this does not solve the problem of spot 
overlap or gross aberration. 

[00012] There have also been a number of other methods for extending the dynamic range of 
the wavefront sensor purely through analysis. For example M. C. Roggeman et al., 
"Algorithm to increase the largest aberration that can be reconstructed from Hartmann sensor 
measurements," Appl. Optics 37(20), pp. 4321^*329 (10 July 1998) proposed a method for 
using image metrics from a separate camera to improve the dynamic range, and Pfund et al., 
"Dynamic range expansion of a Shack-Hartmann sensor by using a modified unwrapping 
algorithm," Opt. Letters 23, pp. 995-997 (1998) proposed a method for analyzing the Shack- 
Hartmann image by use of a modified unwrapping algorithm. While these concepts are 
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useful (in fact, one such method that has particular advantages is disclosed in detail below), 
they do not solve the problem of spot crossover when large changes in the wavefront gradient 
need to be measured. 

[00013] The '917 patent discloses a method for the measurement of the eye that includes, 
among other elements, an adjustable position optical system that provides for an extension of 
the dynamic range of a wavefront sensor for measurement of ocular systems. This method 
provided a means for adjusting a reference sphere through the movement of one optical 
element relative to another optical element, and thus providing a means for limiting the 
dynamic range of the wavefront that is incident upon the sensor. A means for finding the 
appropriate location of this sphere was provided through feedback from the sensor. 
[00014] This system has the advantage of incorporating a very large dynamic range, while 
still providing excellent accuracy. In this instrument, the defocus term was subtracted 
optically so that the residual aberrations were within the dynamic range of the wavefront 
sensor. However, it was applicable primarily to ocular systems that i) permitted injection of a 
small diameter beam, and ii) had well separated aberrations where focus dominated the 
strength of other aberrations. For an arbitrary aspherical optic, these features are not 
necessarily present. 

[00015] Accordingly, it would be advantageous to provide a system and method of 
measuring and mapping three-dimensional surfaces of spherical and aspherical objects. It 
would also be advantageous to provide such a system and method that operates with an 
improved dynamic range. Other and further objects and advantages will appear hereinafter. 
[00016] The present invention is directed to a system and method of measuring and mapping 
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[00017] In one aspect of the invention, a variable position optical system and a dynamic- 
range-limiting aperture are used, similar to that disclosed in the '917 patent, to ensure that a 
wavefront sensor operates always within its dynamic range. It is recognized that for an 
arbitrary optical device with strongly varying surface curvatures, it will not be possible to 
measure the entire element in a single operation (as it is for the ocular system using the 
system disclosed in the '917 patent). However, the inventors have recognized that, with 
different positions of the variable position optical system, it is possible to examine the entire 
surface of the element under test. Thus by systematically varying the position of the variable 
position optical system, it is possible to obtain measurements of the entire surface of the 
element under test, and then construct the measurement of a highly curved or aberrated 
element from these multiple individual measurements. 

[00018] In another aspect of the invention, a system for controlling, directing, and 
modulating light is used to project light onto an element under test, collect the light reflected 
from this element, direct light through a dynamic-range-limiting aperture, and project this 
light onto a wavefront sensor. This system may include optical reformatting optics to 
appropriately reform the light from collimated to converging or diverging as needed for the 
element under test. 

[00019] In another aspect of this invention, a series of wavefront measurements are 
"stitched" together using mathematical methods. Each of these measurements would be 
acquired using a different optical system aspect (such as focus, tilt or other systematically 
introduced reference aberration), in conjunction with a means for limiting the dynamic range 
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of the wavefront incident on a wavefront sensor (beneficially, through a dynamic-range- 
limiting aperture), so that a series of accurate, independent measurements are acquired from 
the element under test. Using the a priori reference information, each individual 
measurement wavefront is corrected appropriately. These measurements may then be 
combined together to produce an overall measurement of the entire surface of the element 
under test using the mathematical methods disclosed herein. 

[00020] Accordingly, in one aspect of the invention, a system for mapping a surface of a 
three-dimensional object, comprises: a projecting optical system adapted to project light onto 
an object; a pre-correction system adapted to compensate a light beam to be projected onto 
the object for aberrations in the object, the pre-correction system being positioned in between 
the projecting optical system and the object; an imaging system adapted to collect light 
scattered by the object; and a wavefront sensor adapted to receive the light collected by the 
imaging system. 

[00021] In another aspect of the invention, a method of mapping a surface of an object, 
comprises: projecting a light beam onto an object; compensating the light beam to be 
projected onto the object for aberrations in the object; collecting light scattered by the object; 
and sensing a wavefront of the collected light scattered by the object. 
[00022] In yet another aspect of the invention, a system for measuring an optical 
characteristic of an optically transmissive object, comprises: a projecting optical system 
which projects light through an optically transmissive object; a correction system adapted to 
at least partially compensate a light beam that has been projected through the object for at 
least one optical property of the object; an imaging system adapted to collect the light that has 
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been projected through the object; and a wavefront sensor adapted to receive the light 
collected by the imaging system. 

[00023] In another aspect of the invention, a method of measuring an optical quality of an 
optically transmissive object, comprises: projecting a light beam through an optically 
transmissive object; at least partially compensating the light beam that has been projected 
through the object for at least one optical property of the object; collecting the light beam that 
has been projected through the object; and sensing a wavefront of the collected light. 
[00024] In another aspect of the invention, a method of mapping a surface of an object, 
comprises: (a) projecting a light beam onto a surface of an object; (b) collecting light 
scattered by a first portion of the surface of the object and rejecting light scattered by a second 
portion of the surface of the object; (c) sensing a wavefront of the collected light returned by 
the portion of the surface of the object; (d) repeating steps (a) through (c) for a plurality of 
different portions of the surface of the object that together span a target area of the surface of 
the object; and (e) stitching together the sensed wavefronts to produce a complete 
measurement of the target area of the surface of the object. 

[00025] In another aspect of the invention, a method of measuring an optically transmissive 
object, comprises: (a) projecting a light beam through a portion of an object; (b) collecting 
light passed through the portion of the object; (c) sensing a wavefront of the collected light 
passed through the portion of the object; (d) repeating steps (a) through (c) for a plurality of 
different portions of the object that together span a target area of the object; and (e) stitching 
together the sensed wavefronts to produce a complete measurement of the target area of the 
object. 
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[00026] In another aspect of the invention, a method of mapping a surface of an object, 
comprises: (a) locating a light source a first distance from an object; (b) projecting a light 
beam from the light source onto a surface of the object; (c) collecting light scattered by the 
surface of the object; (d) sensing a wavefront comprising a difference between a wavefront of 
the collected light and a reference wavefront; (e) changing the distance between the light 
source and the object; (f) repeating steps (b) through (e) to produce N sensed wavefronts; and 
(g) stitching together the N sensed wavefronts to produce a complete measurement of the 
target area of the surface of the object. 

[00027] In another aspect of the invention, a method of measuring an optically transmissive 
object, comprises: (a) locating a light source a first distance from an optically transmissive 
object; (b) projecting a light beam from the light source through the object; (c) collecting light 
projected through the object; (d) sensing a wavefront comprising a difference between a 
wavefront of the collected light and a reference wavefront; (e) changing the distance between 
the light source and the object; (f) repeating steps (b) through (e) N times to produce N sensed 
wavefronts; and (g) stitching together the N sensed wavefronts to produce a complete 
measurement of the target area of the surface of the object. 

[00028] In another aspect of the invention, a point light source for producing a spherical 
wave, comprises: a light source; a diffuser adapted to receive light from the light source; and 
a structure having an aperture adapted to receive and pass therethrough the light from the 
diffuser. 

[00029] In another aspect of the invention, a method of determining when a portion of a light 
wavefront received by a wavefront sensor exceeds the dynamic range of the wavefront sensor, 
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comprises: assigning a group of N pixels of a wavefront sensor to a focal spot; providing a 
first light wavefront to the wavefront sensor under conditions known to be within a dynamic 
range of the wavefront sensor; calculating a reference value, Ok REF , for a second moment of 
the focal spot produced by the first light wavefront within the group of N pixels; providing a 
second light wavefront to the wavefront sensor; calculating a value of the Ok, for a second 
moment of the focal spot produced by the second light wavefront within the group of N 
pixels; and determining that the second light wavefront is within the dynamic range of the 
wavefront sensor within the group of N pixels when: \cr k - \<t a , where t D is a set 
threshold value. 

[00030] In another aspect of the invention, a method of mapping a surface of an object, 
comprises: projecting a light beam onto an object; compensating the light beam to be 
projected onto the object for aberrations in the object; passing light scattered by the object 
through a dynamic-range-limiting aperture; collecting light passed through the dynamic- 
range- limiting aperture; and sensing a wavefront of the collected light. 

[00031] In another aspect of the invention, a method of determining a position of a focal spot 
on a wavefront sensor, comprises: assigning a first group of N pixels of a wavefront sensor to 
a focal spot; providing a light wavefront to the wavefront sensor; measuring an irradiance 
distribution of the light wavefront across the N pixels of the first group; calculating a 
preliminary centroid position of the focal spot as a first moment of the irradiance distribution 
of the light wavefront across the N pixels of the first group; assigning a second group of N 
pixels of the wavefront sensor to the focal spot, where the second group of N pixels is 
approximately centered at the preliminary centroid position; and calculating a location of the 



10 



Docket No. WFS.006CIP 
- PATENT 



focal spot as a first moment of the power of irradiance distribution of the light wavefront 
across the N pixels of the second group. 

[00032] In another aspect of the invention, a method of determining a wavefront of light 
received by a wavefront sensor, comprises: (a) providing a light wavefront to a wavefront 
sensor; (b) assigning pixels of the wavefront sensor to a first plurality of areas-of-interest 
(AOIs); (c) determining a first region of the wavefront sensor where the received light 
wavefront is within a dynamic range of the wavefront sensor for all AOIs within the first 
region; (d) calculating locations for centers of light spots of received light for all AOIs within 
the first region; (e) calculated a fitted wavefront for the received light wavefront over the first 
region; (f) computing a slope of the fitted wavefront at each AOI within the first region; (g) 
projecting locations for centers of light spots of received light for a second region of the 
wavefront sensor larger than the first region, using the slopes of the fitted wavefront within 
the first region; (h) reassigning the pixels of the wavefront sensor to a second plurality of 
areas-of-interest (AOIs) each centered on one of the calculated or projected centers of light 
spots; (i) determining a new first region of the wavefront sensor where the received light 
wavefront is within a dynamic range of the wavefront sensor for all AOIs; and (j) repeating 
steps (d) through (i) until one of: (i) the new first region is no larger than a previous first 
region; and (ii) the new first region spans substantially the entire wavefront sensor. 
[00033] In another aspect of the invention, a method of measuring a focal length (F) of a 
lens, comprises: (a) locating a light source on a first side of a lens, one of the light source and 
the lens being located at a position Zi; (b) locating a wavefront sensor a fixed distance (L) 
from the lens on a second side thereof; (c) projecting a light beam from the light source 
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through the lens; (d) collecting light passed through lens; (e) sensing a wavefront of the 
collected light at the wavefront sensor; (f) measuring a corresponding vergence Pi of the light; 
(g) incrementing i by 1, and moving the position of one of the light source and the lens to a 
new position Zi; (h) repeating steps (c) through (g) to obtain N values of Zi and Pi; and (i) 
applying the N values of Zi and Pi to a least squares fit algorithm to solve for the focal length 
(F). 

BRIEF DESCRIPTION OF THE DRAWINGS 

[00034] FIG. 1 A is a functional block diagram of a system for measuring and mapping three- 
dimensional surfaces of spherical and aspherical objects. 

[00035] FIG. IB is a diagram of an exemplary wavefront sensor that may be used in the 
system of FIG. 1A; 

[00036] FIGs. 2A-B show exemplary positional relationships between certain elements of 
the system of FIG. 1. 

[00037] FIG. 3 shows an exemplary condition where a device under test has a portion of its 
surface that produces reflections that would typically be out of range for a wavefront sensor; 
[00038] FIGs. 4A-E illustrate measurements of portions of a surface of a device under test 
and then stitching together the partial measurements to produce a complete measurement of 
the surface; 

[00039] FIGs. 5 A-D illustrate the effect of out-of-range conditions on the pixels values in a 
group of pixels assigned to a focal spot (Area of Interest). 
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[00040] FIG. 6 is a block diagram description of an algorithm to expand the dynamic range 
by use of iterative wavefront fitting method for tracking Areas of Interest. 
[00041] FIG. 7 shows the measured refractive power (wavefront curvature) as a function of 
stage 1 position. 

[00042] FIG. 8 shows an arrangement for measuring transmissive optical elements using the 
methods of this invention. 

[00043] FIG. 9 shows an arrangement for measuring a transmissive optical element that has 
significant optical power. 

[00044] FIG. 10A shows an arrangement of elements to produce a point source that can be 
used for transmissive or reflective optical measurements. 

[00045] FIG. 10B shows an arrangement of elements to produce a point source of light with 
high power for optical metrology measurements 

[00046] FIG. 1 1 shows the parameters needed for measuring the focal parameters of a lens. 
[00047] FIG. 12 shows example data from the measurement of lens focal properties. 

DETAILED DESCRIPTION 

[00048] FIG. 1 shows a functional block diagram of one embodiment of a system 100 for 
measuring and mapping three-dimensional surfaces of spherical and aspherical objects. 
[00049] The system 100 includes: a projection system 110 comprising a light source 1 12, a 
collimating lens 114, and a polarizing beam splitter 116; an optical imaging system 120 
comprising lens pair 122 and 124 (operating together as a telescope 125), dynamic-range- 
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limiting aperture (RLA) 115, \/4 waveplate 126, and reformatting lens 129; a wavefront 
sensor 130; and a data analysis system 150 including a processor 152. Beneficially, the 
system also includes one or more movable stages - in particular the system 100 includes three 
movable stages 105, 107, and 109 - on each of which one or more components may be 
mounted. The purposes and function of the movable stages will be explained in further detail 
below. 

[00050] Beneficially, in one embodiment, the wavefront sensor 130 is a Shack-Hartmann 
wavefront sensor, an example of which is shown in FIG. IB. The Shack-Hartmann wavefront 
sensor of FIG. IB includes a lenslet array 132 that dissects the incoming wavefront 133 and 
projects focal spots 135 onto a sensor (e.g., a detector array) 134. Although the Shack- 
Hartmann wavefront sensor provides some benefits in certain contexts, as explained in detail 
in various places below, the system 100 of FIG. 1A is not limited to the use of a Shack- 
Hartmann wavefront sensor and, instead, other sensors, such as a shearing interferometer, 
may be used, without departing from the broad scope of various aspects of the invention. 
[00051] It should be noted that the term dynamic-range- limiting aperture (RLA) refers to an 
aperture for limiting the wavefront dynamic range of light that is incident upon the wavefront 
sensor. In most cases this may be an aperture with a prescribed shape that is located relative 
to an imaging lens in the optical system so that it limits the range of angles that may pass 
through the system. The RLA size and shape is adjusted so as not to exceed the range that the 
wavefront sensor can detect without ambiguity. Beneficially, the RLA has a rectangular 
shape (e.g., a square shape) to match the shape of the pixel group assigned to a lenslet, for 
example. However, the RLA may have a round shape, hexagonal or other shape as needed by 
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the particular sensor. The RLA may be designed to pass light to exceed certain strict limits of 
the dynamic range to provide for an optimum combination of optical and algorithmic 
analysis. Such an approach may allow the actual dynamic range of the system to be 
increased. 

[00052] FIG. 1A also shows an exemplary device under test (DUT) 170. The DUT has a 
non-planar, three-dimensional surface 175. In the example shown in FIG. 1A, the surface 175 
is generally spherical in nature, and more particularly, may be a mold for manufacturing a 
contact lens, intraocular lens, aspheric lens or other optical element. It may further be a 
telescope mirror, aspheric imaging lens, lens element, or photolithography lens, or any other 
surface that has specular reflectivity. The system 100 is able to measure and map the surface 
175 to characterize any surface irregularities and aspherical aberrations as described in more 
detail below. 

[00053] FIGs. 2A-B help to explain the positional relationship between the telescope lens 
124, the reformatting lens 129, and the DUT 170. The projection system 1 10 is not shown in 
FIGs. 2A-B for simplicity of illustration. As illustrated in FIG. 2A, for a concave DUT 170 
light passing through the telescope lens 124 is imaged at an image plane position designated 
"13." The reformatting lens 129 is located a distance Si3 from the position of the image plane 
13. The reformatting lens 129 has a focal length "B" and the surface 175 of the DUT 170 has 
a base radius of curvature that is indicated as "R" in FIG. 2A. 

[00054] In this case, given the base radius of curvature R of the surface 175 of the DUT 170 
and the focal length f3 of the reformatting lens 129, one can calculate the position Si3 to be: 



15 



Docket No. WFS.006CIP 
-- PATENT -- 



1 1 (i) 



Si, f, (/,+*) 



[00055] FIG. 2B illustrates the case of a convex DUT 170. In this case the location of the 
image plane 13 is adjusted so that it matches the virtual image of the DUT 170 created by the 
reformatting lens 129 (L3). In this case the radius R is negative and the resulting image 
distance Si3 is therefore also negative. By adjusting the position of both reformatting lens 
129 and the DUT 170 (i.e., by changing the position of the movable stage 107), the image or 
virtual image plane of reformatting lens 129 can still be matched to the object plane (labeled 
13) of telescope lens 124 (LI) and hence imaging can be maintained throughout the system 
from the surface of the DUT 170 onto the wavefront sensor 130. 

[00056] Using either the arrangement of FIG. 2 A or FIG. 2B, highly curved DUTs 170, 
either convex or concave in nature, can be tested. The choice of reformatting lens 129 is 
determined by the radius of curvature of the DUT 170s to be measured. The field of view is 
determined by the magnification M2 = -fl/f2 (where fl is the focal length of lens LI and f2 is 
the focal length of lens L2), and by the size of the wavefront sensor input aperture. 
[00057] Since the device under test may have a wide range of local radii of curvature, it is 
possible that the light that is collected by reformatting lens 129, and transmitted through 
telescope lens 124 (LI in FIGs. 2A-B) to the RLA 115 arrives at a significant angle. If this 
light were to be allowed to propagate through the rest of the system, it would result in a 
wavefront that would exceed the dynamic range of the wavefront sensor 130. The dynamic- 
range-limiting aperture (RLA) 115 clips any such rays, thereby insuring that the wavefront 
sensor 130 never receives light that exceeds its measurement capabilities. However, this 
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means that there may be portions (even significant portions) of the DUT 170 that cannot be 
measured, and hence no light would be returned through the system to the wavefront sensor 
130 for this condition. 

[00058] FIG. 3 shows a DUT 170 that has a portion of the surface that would exceed the 
angular dynamic range of the wavefront sensor 130. A ray that is incident on this portion of 
the wavefront sensor 130 is shown. This ray originates from the injection system, is relayed 
onto the DUT 170, and the collected light is then collected by reformatting lens 129 and 
directed back toward the wavefront sensor 130. However, because this ray has a greater angle 
than the wavefront sensor 130 is capable of detecting, it arrives at the RLA 115 off axis. The 
size of the RLA 1 1 5 is chosen so that only in-range rays may be admitted. Thus the ray 
shown in FIG. 3 is clipped, and does not pass through the RLA 1 1 5 to thereby cause 
inaccurate measurements. The particular ray shown in Figure 3 is from a region that is more 
highly curved than the rest of the DUT 170 (from the inner portion, in this example). 
However, rays from the outer region (in this example) are well matched to the optical system, 
and thus would transit the entire system and would therefore be detected by the wavefront 
sensor 130. All such rays corresponding to the outer region (in this example) would pass 
through the RLA 115 and thus be detected. Thus, in this configuration, the wavefront sensor 
130 would receive light from a ring corresponding to the outer region of the DUT 170. No 
light would be collected from the middle (except for one point exactly on axis). All of the 
light collected in this ring could be used to make an accurate measurement of the DUT 170 in 
this region. 

[00059] This basic concept of using a dynamic-range-limiting aperture was used to good 
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effect in the '917 patent to assure that all rays that were collected by (in this case) the 
ophthalmic wavefront sensor system for measuring the eye were always in range. The 
instrument described in the '917 patent, however, was designed to acquire all of the 
information in a single frame, after adjustment of movable parts of the optical system to 
optimize the dynamic range. The RLA was used, in that case, primarily to ensure that only 
"in-range" light was incident on the sensor during the adjustment procedure, and thus 
facilitate rapid and accurate adjustment. 

[00060] Now, this concept is extended further. By adjusting the position of the moving stage 
105 in FIG. 1 A, the spacing of telescope lenses 122 and 124 may be changed. This changes 
the focus condition for light transmitted through the optics and also received by the wavefront 
sensor 130. Adjustment of the distance between the principle planes of these lenses to a 
distance equal to the sum of their effective focal lengths results in a planar wavefront 
transmitted through telescope lenses 122 and 124. Reformatting lens 129 reformats this to a 
spherical wavefront that converges (in the case of convex DUTs 170) or diverges (for 
concave DUTs 170) to match the appropriate radius of curvature of the DUT 170. If the 
distance between the telescope lenses 122 and 124 is set to a different distance (by adjusting 
movable stage 105), then the transmitted wavefront is no longer planar, but may be 
converging or diverging. Reformatting lens 129 then reformats this to either a slightly longer, 
or a slightly shorter (depending upon the type of DUT 170) radius of curvature wavefront that 
would be incident upon the DUT 170. If this incident wavefront approximately matches the 
radius of curvature of the DUT 170, then the incident light will be reflected back such that it 
can be collected by the optical system, pass through the RLA 115, and be measured by the 



18 



Docket No. WFS.006CIP 
- PATENT - 

wavefront sensor 130. If it does not match the radius of curvature of the DUT 170 well, then 
it will be clipped by the RLA 115. Thus to measure the entire surface of the DUT 1 70, the 
instrument must have sufficient range of travel of the moving stage 105 to allow adjustment 
of the relative positions of telescope lenses 122 (L2) and 124 (LI) over an appropriate range. 
[00061] Accordingly, the telescope lenses 122 and 124, and where necessary the 
reformatting lens 129, function as a correction, or pre-correction, system for that compensates 
the light beam for one or more aberrations in the DUT 170. Alternatively, instead of one or 
more lenses being mounted on a movable stage, a variable focal length lens may be 
employed. The focal length of such a lens may be controlled by the processor 152. 
[00062] At each relative position of telescope lenses 122 and 124, there is a particular radius 
of curvature that can be measured on the DUT 170. Since the DUT 170 may have a wide 
variation in radii, at any one time it will not be possible to measure the entire DUT 170 
(unlike the eye, where the '917 patent instrument was designed to have sufficient dynamic 
range (after adjustment) to measure the entire eye in a single measurement). But a series of 
measurements, each corresponding to a different position of lens 122 relative to lens 124, 
would result in a measurement of different regions of the DUT 170, each corresponding to a 
different radius of curvature. Since the separation between telescope lenses 122 and 124, and 
the focal lengths of these lenses, is known from the instrument design, the relative defocus 
can be determined through the position of the stage 105. Thus the series of measurements, 
each at a known incident radius of curvature can be obtained. 

[00063] Given such a sequence of measurements, it is also possible to take advantage of the 
known incident radius of curvature for each individual measurement, and thus produce the 
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measurement over the entire DUT 170. The sequence of measurements can be pieced 
together, or "stitched" together, to form a single map of the three-dimensional surface of the 
device under test. This is different from what is disclosed in Neal et al. U.S. Patent 6,184,974, 
"Apparatus and Method for Evaluating a Target Larger Than a Measuring Aperture ... in 
that, instead of stitching together small portions of the measurement of the DUT by multiple 
sub-aperture measurements (area stitching), the stitching occurs through stitching together 
portions of the measurement that would otherwise be out of range for the wavefront sensor. 
[00064] The geometry proposed in FIGs. 1-3 for measuring highly curved optical elements is 
not the only geometry where this technique is useful. Modern contact lenses, intraocular 
lenses, and highly aspheric optical elements also pose difficulties in measurement due to 
rapidly changing curvature over part of the optical element. These elements may be better 
measured in transmission, as shown in FIGs. 8 and 9. 

[00065] FIG. 8 illustrates an optical element as the device under test (DUT), such as a 
contact lens 842, which has portions that would normally exceed the dynamic range of the 
sensor 830. FIG. 8 shows measurement of a contact lens or other transmissive optical 
element that has low net power. A collimated beam is passed through the DUT 842 and 
analyzed by the sensor. Either the source 862 or the moving stage 805 can be moved to 
acquire data over the various zones as needed. In this case the lens has little net optical 
power, or at least falls within the range of optical power that may be corrected by moving the 
movable stage 805. Light from a source 862 is collimated using lens 860 and transmitted 
through the DUT 842 which may be held in a test cell or mount 840. This light is collected 
by optical system 825, which contains parts mounted on a moving stage 805. The wavefront 
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sensor 830, dynamic-range-limiting aperture 815 and one of the imaging lenses 816 are 
mounted on this movable stage 805. Advantageously, the wavefront sensor 830 is positioned 
so that the distance from the wavefront sensor 830 to the lens 816 is equal to the focal length 
of lens 816. Meanwhile, the DUT 842 is positioned at a distance from lens 818 equal to the 
focal length of lens 818. During the measurement process the dynamic-range-limiting 
aperture 815 clips any light that would otherwise exceed the dynamic range of the wavefront 
sensor 830. The moving stage 805 position is varied and data is acquired and saved by the 
processor 852. This data may be stitched together using the methods set out below, to 
reconstruct a measurement of the entire DUT 842. Instead of moving the sensor 830, lens 
816, and RLA 815, it may be advantageous to move the source 862 relative to the DUT 842 
instead. This will serve a similar purpose. 

[00066] In a similar manner, an optical element that has significant optical power can be 
measured by omitting the collimating lens 860, and using the DUT itself to partially collimate 
the light from the source. Such a case is shown in FIG. 9. FIG. 9 shows measurement of an 
intraocular lens (IOL) or other transmissive optical element that has significant focusing 
power. A collimated beam is passed through the DUT 942 and analyzed by the sensor. 
Either the source 962 or the moving stage 905 can be moved to acquire data over the various 
zones as needed. In this case the DUT 942 collimates the light from the source 962. This 
light is then collected with imaging system 925 onto sensor 930. In this case there are two 
different ways for the multiple measurements to be stitched together can be obtained. Either 
the position of the movable stage 905 may be varied, as described previously, or the position 
of the source 962 may be varied using movable stage 908. Since the focal power (effective 
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focal length) of the DUT 942 may not be known, it may be necessary to include an option for 
moving the movable stage 908 in any case. Then the movable stage 905 may be omitted 
since it serves no unique function. 

[00067] Now, embodiments of a process for stitching together individual measurements will 
be described. FIGs. 4A-D depicts a sequence of measurements for a particular type of DUT 
170 using the system of FIG. 1 A. In this case the DUT 170 is a mold for a bifocal contact 
lens. The mold has regions with widely different radii of curvature with extremely sharp 
transitions between regions. As can be seen in FIG. 4E, there are six basic annular zones, 
with zones 1, 3 and 5 having one radius of curvature, zones 2 and 4 a different radius of 
curvature, and zone 6 having yet a different radius. It is important to be able to measure the 
different radii of curvature accurately, e.g. an absolute accuracy of better than 0.005 mm for a 
DUT 170 with a 7 mm base curve. For the data represented in FIGs. 4A-D, the radii of 
curvature of the different regions is sufficiently different that it would not be possible to 
design a wavefront sensor 130 (e.g. a Shack-Hartmann wavefront sensor) with sufficient 
dynamic range that could also achieve the required accuracy. 
[00068] The following techniques are applied to measure the DUT 170. 
[00069] First, the DUT 170 is positioned in the metrology instrument 100 at the appropriate 
image plane. This plane may be determined by either optical or mechanical means. 
[00070] To determine the appropriate position optically, the DUT 1 70 is initially positioned 
with the surface approximately one focal length from reformatting lens 129. The telescope 
lenses 122 and 124 are separated to be spaced apart by the sum of their focal lengths, so that 
they will provide a collimated (flat) transmitted wavefront when the DUT 170 is positioned 
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exactly at the image plane. Then, the position of the DUT 170 is adjusted (by means of 
movable stage 109) until the wavefront as measured by the wavefront sensor 130 is flat, 
which means that the DUT 170 is located at the so-called "cat's eye" position. This is a well- 
known condition optically where the light is focused on the surface of the DUT 170 at a 
single point. It is not sensitive to alignment or position of the DUT 170 other than the 
relative distance from reformatting lens 129. This position is recorded using a micrometer, 
encoder or other position gauge on the relative position between the DUT 170 and 
reformatting lens 129. 

[00071] Next, the DUT 170 is moved (again, by movable stage 109) by approximately the 
radius of curvature of the central region of the DUT 170. The position of the DUT 170 is 
then finely adjusted to once again provide a flat wavefront over at least some portion of the 
surface (usually the central region) as determined by the wavefront sensor 130. This is often 
called the "confocal position." The difference between the "cat's eye" position, and the 
confocal position is the base radius of curvature of the DUT 170 (Rb). 

[00072] If the DUT 1 70 is arranged so that it may be registered accurately (at least to within 
the desired measurement accuracy), then the preceding steps (moving to the "cat's eye" 
position and then the confocal position) may not be necessary once the instrument is 
calibrated. The calibration and positioning accuracy are sufficient to achieve accurate 
measurement and the base radius may be determined only from the stage position readout at 
the collimation position. 

[00073] Next the spacing between the imaging lenses 125 is varied systematically so as to 
acquire a sequence of measurements. In the example that produced the chart of FIG. 4E, 40 
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separate measurements were acquired. However, hundreds or even thousands of 
measurements may be used to increase the signal-to-noise ratio. Each measurement is taken 
with a different spacing between telescope lenses 122 and 124 and thus with a different 
reference defocus. The spacing is accurately recorded using a position detection system (e.g. 
an encoder) on the mechanical stage that is used for positioning. For a wavefront sensor 130 
that measures wavefront slope (such as Shack-Hartmann, Moire deflectometry, shearing 
interferometer), the position of the stage 105 can be used to modify the acquired wavefront 
data to include the appropriate defocus. 

[00074] The surface shape, S(x, y), is determined partly by the defocus introduced through 
the variable separation of telescope lenses 122 and 124 and reformatting lens 129 and hence 
the stage position f(x,y), and partly from the measured wavefront of the sensor w(x,y): 



(2) 



[00075] The resulting slope at each point across the aperture is: 



df(x,y) 
dx 
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Where the defocus wavefront is given by f(x,y): 
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and so 
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[00076] Each of these measurements may not have a complete characterization over the 
surface, but the collection together should span the entire space. If the measurements are 
from a slope type wavefront sensor (as in the example above) then the slope distributions 
must be reconstructed to retrieve the wavefront distribution. 

[00077] FIGs. 4A-D illustrate a sequence of these measurements. FIGs. 4A, B, and C show 
three different individual measurements. Note that there is a significant portion of each 
image where no data was acquired. However, these regions overlap in such a way that the 
entire image can be retrieved. 

[00078] The images are stitched together using Eqns. 2 - 6 to include the appropriate 
defocus term. The stitching can be done through a number of methods, which may include 
spot tracking, averaging, or other means. 

[00079] As an example, averaging can be applied to defocus corrected measurements. In the 
averaging method for stitching, one takes advantage of the fact that Eqns. 2-6 apply the 
correct defocus term to each of the wavefront sensor measurements so that each measurement 
represents a complete, accurate measurement of the surface, even if only a portion of the 
surface has available data. Thus for each location (x,y) on the surface, the surface height at a 
given position is determined independently by each measurement. Thus, an accurate measure 
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of the surface height can be determined simply by averaging each surface height from each 
valid measurement. However, for strongly varying DUT 170s, there are regions where no 
valid data is available (due to the RLA 115). These points are not included in the average, so 
that for each location (x,y), a different number of points may be included in the average, 
depending upon the number of valid overlapping measurements. As the defocus is varied, 
this has the effect of "filling in" the missing portions of the data so that the whole surface is 
included. 

[00080] For a slope type wavefront sensor, the wavefront is determined by integration from 
wavefront gradient measurements. This integration involves an arbitrary constant of 
integration that cannot be specified. So, for unconnected regions it is impossible to determine 
the actual surface position, and only relative surface height within a connected region can be 
obtained. However, the defocus correction of Eqns. 5 and 6 applies to the slope. Therefore 
if, instead of averaging the wavefronts, the wavefront slopes are averaged, then it is possible 
to perform the wavefront reconstruction after the stitching process has completed. This will 
lead to a large connected region and hence improved performance of the algorithms. 
Mathematically, for an AOI "k" and measurement "m," the slope at a point 0 k x (x,y) can be 
constructed from a series of measurements 0 kjm x (x,y) given a validity function V k , m (x,y), 
similarly for 9k, m y (x,y). 



The validity function describes where data is valid and may easily be obtained by evaluating 
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the irradiance of the light at each measurement point. Among other methods, the validity 
function can be constructed from irradiance measurements I k , m (x,y). 



where t is a predetermined "irradiance threshold." 

[00081] In the embodiment where the wavefront sensor 130 consists of a Shack-Hartmann 
wavefront sensor that comprises a lenslet array disposed in front of a detector array, there are 
a number additional ways to construct the validity function in an optimal fashion. The total 
light that is collected by each lens in the lenslet array is readily determined by adding up all 
the pixel values associated with that lenslet. This can be used as an estimate of the irradiance 
for this portion of the measurement Ik(x,y). Eqn. 8 can be used to construct the validity 
function through the use of the RLA 1 15 to clip any light that would otherwise be out of 
range. This validity information may be obtained in other ways from different measurement 
instruments. However, it is usually easy to determine a region where valid data exists by 
analyzing the images obtained by the respective sensor. 

[00082] However, additional information may also be obtained from the Shack-Hartmann 
image. Each focal spot covers a number of pixels. This is necessary in order to get optimum 
accuracy from this type of sensor. Using a Shack-Hartmann wavefront sensor and sensing 
method, an area-of-interest (AOI) comprising N pixels is assigned to each lenslet. Neal et al., 
"Shack-Hartmann wavefront sensor precision and accuracy," SPIE 4779, (2002), showed that 
that optimal signal to noise was obtained when about 50 pixels were covered by each focal 
spot in the measurement process. This is about a 7 X 7 group of pixels, in general. The 
position of the focal spot is determined through a 1 st moment algorithm (sometimes called a 
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centroid algorithm): 



x = 



i 



(9) 
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where the summation is performed over all N pixels in each area-of-interest (AOI) assigned to 
each lenslet (similarly for the y position). With about N=50 illuminated pixels, the focal spot 
position can be determined to about 1/5 0 th of a pixel, which is the boundary where other 
systematic effects start to affect the accuracy. With this distribution of pixels, however, it is 
also possible to estimate the size of the focal spot. This can be done through the 2 nd moment 
algorithm: 



[00083] FIGs. 5 A-D illustrate several exemplary conditions of illuminated pixels in a single 
AOL 

[00084] FIG. 5 A shows a typical AOI with normal in-range light. The focal spot size comes 
into play because it provides information about the changes in the lenslet focal spot. The 
lenslet focal spot can be modified by various effects in the optical system which it would be 
advantageous to detect. 

[00085] FIG. 5B illustrates a slightly out-of range condition. Some light from an adjacent 
lenslet is "leaking" into the edge of the pixels under the lenslet of interest. This biases the 
data and reduces the accuracy. In this condition the 2 nd moment spot size changes 
dramatically (by factors of about 4-10 or more) since the positional weighting (x 2 ) gives a 
strong weight to a small amount of signal on either edge of the area of interest. A very coarse 
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threshold (i.e. 2-3 times the average a (averaged over all k)) is sufficient to detect this 
condition and add these AOIs to the validity function, i.e. 
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where the reference second moments, o k , are recorded in a calibration step where the 
collected wavefront was ensured to be completely in range. Usually this is the result of 
calibration with a perfect sphere or other ideal optical element. 

[00086] FIG. 5C illustrates a diffracted focal spot. Since the operation of the system as 
disclosed above necessitates scanning the defocus position, it is evident that the focal spots 
will move though the AOIs, entering from one edge and leaving from the other. This will 
perforce cause the focal spot to transit the RLA 115. When the focal spot is partly on the 
edge of the RLA 115, this clipping will cause diffraction. This diffraction will cause the focal 
spot to grow until it is completely clipped. Using Equations 10 and 1 1, it is possible to 
identify these conditions and construct the appropriate validity function. 
[00087] FIG. 5D illustrates a case where the DUT 170 has a complex geometry. For the 
complex DUT 170 described above, there are regions with a sharp transition in curvature. 
This means that there are cusps or sharp changes in slope on the surface of the DUT 170. For 
those lenslets that cover such a sharp transition, the light will be reflected in two different 
directions simultaneously. This will result in greatly distorted focal spots. These distorted 
focal spots will be bigger than a non-distorted focal spot, and hence can be readily detected 
with Eqns. 10 and 11. 
[00088] Focal spot tracking. 
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[00089] For a Shack-Hartmann wavefront sensor, the assignment of pixels to lenslets is 
required to assure that the wavefront can be determined without ambiguity. Thus a unique 
mapping from a particular group of pixels to a particular lenslet must be maintained. Several 
different techniques for determining this mapping have been used (e.g., T. Bruno et al. 
"Applying Hartmann wavefront sensing technology to precision optical testing of the Hubble 
Space Telescope correctors," SPIE 1920, pp. 328-336 (1993)). In practice, the most common 
is to predetermine this mapping from an initially flat wavefront that is incident upon the 
sensor. In this case the dynamic range of the sensor is given by the location where the edge of 
the focal spot is located at a position that is on the boundary of the pixels assigned to that 
lenslet. For a sensor designed without intervening optics between the lenslet array and the 
detector array (e.g., following the methods of U.S. Patent 6,130,419 "Fixed mount wavefront 
sensor"), the size of the group of pixels is the same size as the lenslet array. Thus the angular 
dynamic range is given by 

e _ m 41=£I± (12) 

[00090] In practice, the use of thresholds means that this maximum angular dynamic range 
estimate is usually slightly conservative, that is, the threshold focal spot is actually slightly 
smaller than fX/d in radius, so the second term in Eqn. 12 is slightly less, resulting in a 
slightly larger dynamic range. For a 12% threshold, the threshold spot size is approximately 
3X/4d. 

[00091] However, it should be noted that the group of pixels assigned to the focal spot is 
purely an arbitrary convention. Thus, using a priori, internal, or other information to increase 
the dynamic range of the sensor can adjust this assignment. 
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[00092] Further, the group of pixels that are assigned to a particular lenslet can be 
determined through an algorithm that assigns the pixel group to the lenslet array, and then 
determines the centroid (or spot location) from the assigned group of pixels. To this end, the 
following algorithms which may be used to assign the pixel groups, or Areas of Interest 
(AOIs), to the lenslets. 
[00093] One frame tracking. 

[00094] If the focal spots are arranged so that there is sufficient local tilt in the wavefront 
(and either no RLA is employed, or the size of the RLA is large) such that two partial focal 
spots are located on the edges of the pixel AOI, then there will be a significant error in 
estimating the spot location using Eqn. 9. However, if the local tilt is less than d/(2f), that is, 
if the centroid as determined by Eqn. 9 is in the correct half of the AOI, then an additional 
calculation step will result in increased dynamic range. This calculation consists of moving 
the pixel group (AOI) to be centered on the centroid position determined by the first step, and 
then calculating again the centroid using Eqn. 9, but with the new group of pixels. Unless the 
focal spots on each edge are precisely balanced, this will result in almost a 50% increase in 
the dynamic range. To assure that the focal spots are in the right arrangement initially, the 
RLA size can be adjusted to clip light that has an incident angle exceeding d/(2f). As an 
example consider a lenslet with 2 mm focal length, and .072 mm diameter. In this case the 
nominal dynamic range (for 635 nm light) is 9.2 mr, whereas one frame tracking would allow 
18 mr of angular dynamic range, almost a factor of 2 improvement. 
[00095] Iterative wavefront fitting. 

[00096] For some types of optical system the wavefront may have some region where it is 
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bounded and in range in some locations, but exceeds the nominal dynamic range of the 
system according to the definition given in Eqn. 12. If the wavefront is smoothly varying, 
then it is possible to extrapolate the location of AOIs outside the region of in-range data. It is 
typical that the in-range portion of the wavefront is in the center of the sensor area, although it 
may be in a different location, as long as this location can be identified using a priori, 2 nd 
moment, or other methods. 

[00097] This method is outlined in FIG. 6. Starting with the region of in-range data, a 
wavefront surface is fit to the data over only the region that is confirmed to be in-range. This 
can be checked using the 2 nd moment algorithm to verify that only a single focal spot is 
contained in each AOI in a particular region. The wavefront fitting is performed using 
Zemike polynomials, Taylor monomials, or other polynomial or continuous functions as 
appropriate to the data under test. A priori information can be used to limit the region of 
connected data if it is available. The fitted wavefront can then be used to project the location 
of the center of an extrapolated set of AOIs by computing the slope at center of each lenslet, 
and projecting the location of the centroid associated with the lenslet by assumption that the 
wavefront is continuous and varies smoothly over at least a small region. The AOI positions, 
that is, the particular pixels assigned to a lenslet, are changed to reflect this expected focal 
spot location, and then Eqn. 9 is used to compute a new centroid. The 2 nd moment algorithm 
(Eqn. 10) can be compared to a threshold to determine if the algorithm is successful in 
calculating in-range focal spots. The resulting centroid pattern is then used to compute a new 
set of wayefront slopes, which are then used in the computation of a new wavefront 
polynomial fit. The process can be repeated until either all of the data has been analyzed, or 
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regions are identified where the data is so rapidly varying that no successful wavefront can be 
computed. 

[00098] Calibration steps. 

[00099] The system 100 of FIG. 1A can be calibrated internally through the use of a few 
known parameters. Since the optics are mounted on movable stages with precise motion 
control, the position of various groups of optics may be varied systematically while recording 
the outputs of the wavefront sensor. The relative position of various groups of optics can be 
precisely determined through the use of common linear travel stages with stepper (or servo) 
motor drives and position encoders. In practice it is easy to obtain accurate positioning on the 
order of 1 micrometer. Variation of the optics position, in various groupings, allows the 
determination of most of the optical parameters of the system. 

[000100] Parameters that need to be determined include, but may not be limited to, the 
focal lengths of the various lenses fi, f 2 , and f 3 , the offset positions of the various stages (that 
is the absolute position of the center of travel of each movable stage), the system 
magnification, and the conversion between stage position and power of equivalent defocus in 
diopters when the stage positions are varied. 

[000101] FIG. 7 illustrates an example of a determination of a cal factor, Ci, for the stage 
105. FIG. 7 plots measured refractive power (wavefront curvature) as a function of stage 1 
position (stage 105). This curve is linear, with a slope 0.3823 diopter/mm. The stage cal 
factor of 2.6155 mm/diopter allows the stage position to be used to correct the optical 
wavefront by a known systematic curvature. For any stage position Zi, the wavefront 
curvature 1/Ri can be determined by just multiplying the stage position by the cal factor: 
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1/Ri=ZiCi. So, Ci is used to determine the correct wavefront curvature to add to the 
measured wavefront during the focal stitching process as described by Eqn. 4. 
[000102] Other factors may determined by systematic variation of the other stages as 
appropriate. 

[000103] One of the key parameters that needs to be determined for proper analysis and 
calibration of the system is the focal length of the various lenses. While there are a number of 
different techniques for determining the focal length of a lens, many of these techniques 
require a subjective evaluation. Other techniques rely on mechanical measurements of the 
lens surface, which make it difficult to account for index of refraction variations or 
aberrations. Even image-based techniques rely on determination of best focus to determine 
the focal properties of the lens. For lenses with variations this leads to inaccuracies in 
determining the focal length. 

[000104] The wavefront sensor described previously can measure the radius of curvature of 
an incident wave. This can be used to determine the focal properties of a lens in the 
following manner, as explained with respect to FIG. 1 1 . FIG. 1 1 illustrates the geometry for 
measuring lens focal lengths and other properties of lenses using a wavefront sensor. In FIG. 
1 1, a lens to be tested 1 140 is arranged between a point source 1110 and a wavefront sensor 
1 130. The lens 1 140 collects the light from the source 1110 and directs it toward the sensor 
1 130. The sensor 1 130 measures the wavefront of the light incident upon it, including the 
wavefront curvature. The focal length is a paraxial property of the lens 1 140. That is, the 
focal length is the distance that the lens 1 140 will focus the near on- axis (paraxial) rays from 
a collimated source. The focal length is related to the object and image positions by the thin- 
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lens formula: 

7-tT 

[000105] With the wavefront sensor 1 130 arranged as shown in FIG. 1 1, the parameters S 0 
and Si can be replaced with the appropriate quantities from the optical geometry: 

!-— ! +^ d4) 

/ f + Z-Z 0 L + R 

which can be rewritten in terms of the measured vergence at the wavefront sensor P=J/R: 



1 = Z-Zq 
R f 2 +(f-L)(Z-Z 0 ) 



where (Z - Zo) is the relative position of the source (Z) with respect to its position at 
collimation (Zo),/is the lens focal length of lens 1 140, and L is the distance between the lens 
1 140 and wavefront sensor 1 130. This formula expresses the relationship between the 
position of a point source 1110 and the measured power (1/R) of the wavefront incident on 
the wavefront sensor 1 130. If the source 1 1 10 is located on a moving stage that allows for 
precise variation of its position in a known fashion, e.g. through the use of a motorized stage 
with encoder position feedback (such as in the system 100 of FIG. 1 A), then a data set of 
measured power P as a function of source positions can be obtained. This data can be 
analyzed to determine the parameters f, L, and Z 0 that provide the best fit to Eqn. 15. A least 
squares, bisection, iterative or other numerical method can be employed. 
[000106] FIG. 12 illustrates measurement of the focal power (P=l/R) as a function of stage 
position over a range of positions. Analysis of this curve allows the determination of the lens 
focal length and other properties. FIG. 12 presents an example, which shows the measured 
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power (in diopters) as a function of stage position over a +/-0.12 m range. Analysis of this 
data yields f=301.52 mm, L=405.65 mm, and Z 0 =0.21 mm. 

[000107] The analysis of paired object/image location data is well known in the prior art. 
For example, Pernick et al., "Least-squares technique for determining principal plane location 
and focal length," Appl. Optics 26(15), pp. 2938-2939 (1 Aug 1987) showed that lens focal 
lengths could be determined through measurements of paired data and a least squares fitting 
routine. 

[000108] However, in the prior art the image position is determined by examining the image 
and finding the position that results in best focus. This is difficult for an aberrated lens, or for 
one with a long focal length. In the case of an aberrated lens there is no single clear position 
for best focus, and for a long focal length lens there is a long region where the image appears 
to have similar quality because of the large depth of focus for such a lens. In contrast to the 
prior art, in this method a wavefront sensor (Shack-Hartmann, Moire deflectometry, 
holographic grating, shearing interferometer, etc.) is employed. This wavefront sensor can be 
used to describe the incident wavefront through a polynomial decomposition, such as Taylor 
polynomials, Zernike polynomials or other orthogonal or non-orthogonal basis. The 
advantage of using an orthogonal basis (such as Zernike polynomials) is that the paraxial 
effects (such as the paraxial focal power) can be separated from other phenomenon or 
aberrations. Thus by analyzing only the second order polynomial terms, one can analyze the 
paraxial focal power term separately from the higher order terms, such as spherical 
aberrations. This effectively allows the analysis of only the paraxial properties, even in the 
presence of other aberrations. Also, the wavefront determination can be quite accurate, even 
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for extremely large radius of curvature wavefronts. For example, Shack-Hartmann wavefront 
sensors have been shown to be capable of accurately measuring wavefronts having radii of 
curvature >1000 m, over a 5 mm aperture. This accuracy allows the measurement of even 
long focal length lenses that would be difficult using the past. 

[000109] Although the specific example discussed above pertains to a case where the 
absolute position of the device under test (e.g., lens) remains constant while the light source is 
moved by the movable stage to obtain N data samples for the least-squares-mean algorithm, it 
is also possible to move the device under test (e.g., lens) by a movable stage while the light 
source remains fixed, to obtain the N data samples. 

[000110] One common element that is needed for measuring parts using these methods is a 
point source that is used to provide an accurate spherical wave. Such a point source may be 
constructed using a variety of different methods. In the prior art a number of methods have 
been used for such a purpose. These include focusing light through a small pinhole, the use 
of a small optical fiber, or LED. However, each of these methods has serious disadvantages 
that may preclude its use in some situations. Focusing light through a pinhole requires 
accurate and difficult alignment. The numerical aperture from an optical fiber is limited to 
relatively low values (<0.1 1), and an LED may have internal structure that complicate the 
optical wavefront so that it does not present a perfect spherical wave. 
[000111] FIG. 10A illustrates a device for producing a nearly spherical wavefront that 
overcomes these advantages. This device uses a standard LED, which is very low cost and 
yet may provide a significant amount of light power. The LED consists of emitter 1010, with 
contact wire 1030, that are protected inside an epoxy well structure 1020. If only the LED is 
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used as a source, then the contact wire and well structure cause three-dimensional patterns of 
the emitted light waveform that complicate its use as a spherical wave source. By including 
the additional elements of a diffuser 1040 and pinhole structure 1050 with aperture 1060 
these disadvantages may be overcome. The diffuser 1040 may be made of a thin piece of 
Teflon or other appropriate material. 

[000112] Beneficially, when a Shack-Hartmann wavefront sensor is employed, the aperture 
1060 does not need to be extremely small, since it will end up being imaged onto the detector 
array. The size of the effective focal spot may be arranged to grow by only a small amount. 
However, it may be advantageous to arrange a slightly larger spot in any case, since covering 
more pixels will lead to better measurement accuracy (as shown for example by D. R. Neal et 
al., "Shack-Hartmann wavefront sensor precision and accuracy," SPIE 4779, 2002). 
[000113] The disadvantage to this configuration is that the efficiency is very low. The 
aperture 1060 passes only a small amount of light, and this if further reduced by the diffuser 
1040. 

[0001 14] The light efficiency may be improved by using the configuration shown in FIG. 
10B. In this case, a tapered fiber 1070 (or alternatively a tapered fiber bundle) is introduced 
between the LED and the diffuser/pinhole elements 1040/1060. The tapered fiber 1070 
collects the light from the LED, concentrates it, and provides it to the diffuser 1060 and 
pinhole 1070. 

[000115] While preferred embodiments are disclosed herein, many variations are possible 
which remain within the concept and scope of the invention. Such variations would become 
clear to one of ordinary skill in the art after inspection of the specification, drawings and 



38 



Docket No. WFS.006CIP 
- PATENT -- 

claims herein. The invention therefore is not to be restricted except within the spirit and 
scope of the appended claims. 
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